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I.  Summary 

The  research  performed  under  the  contract,  during  the  period  1  May  1084  through 
15  December  1084,  can  be  divided  into  three  major  topics;  modeling  Lg  propagation 
across  regional  transition  zones,  the  use  of  long  period  body  waves  to  constrain  focal 
mechanism  orientation,  and  WKBJ  modeling  of  S  and  SS  to  investigate  the  nature  of 
the  boundary  between  a  tectonic  and  sheld  region. 

In  Section  II,  results  of  modeling  Lg  waves  crossing  continental  margins  by  combin¬ 
ing  analytic  with  finite  element  calculations  are  presented.  The  accuracy  of  the  formula¬ 
tion  is  tested  by  using  a  layer-over-halfspace  model  for  which  a  solution  can  be  indepen¬ 
dently  determined  using  conventional  mode  sum  calculations.  This  test  case  is  also  used 
to  illustrate  the  efficiency  of  an  absorbing  boundary  condition  designed  to  reduce  com¬ 
putation  time.  In  the  modeling  of  transition  regions,  the  decrease  in  crustal  thickness 
from  continents  to  oceans  focuses  seismic  energy.  This  results  in  surface  amplitudes  that 
are  50%  larger  in  the  oceanic  crust  near  the  end  of  the  transition  than  those  in  the  adja¬ 
cent  continent.  However,  the  amplitude  of  the  Lg  waves  rapidly  attenuates  as  they  pro¬ 
pagate  to  larger  distances  in  the  oceanic  crust.  Conversely,  amplitudes  are  reduced  when 
the  Lg  wavetrain  passes  through  the  transition  from  ocean  to  continent  because  energy 
disperses  into  the  thickening  crustal  layer. 

In  Section  m,  a  body-wave  analysis  of  the  1966  Gisborne,  New  Zealand,  earthq^ke 
is  given.  The  ML  =*  6.2  event  occured  along  the  east  coast  of  the  North  Island.  Modeling 
of  P  and  S  body  waves  shows  that  the  focal  mechanism  of  this  event  is  consistent  with 
the  northwestward  thrusting  of  the  Pacific  Plate  beneath  North  Island  (^=-=249*  ,6=25* 
and  X~131* ).  The  focal  depth  is  constrained  to  18  km,  significantly  less  than  the  values 
of  25  to  30  km  computed  from  local  network  data.  Estimates  of  scalar  moment,  source 
duration  and  stress-drop  for  the  event  are  4x10*  dyne-cm,  2  to  3  seconds,  and  20  to  120 
bars  respectively.  Cross-correlation  errors  of  synthetic  to  observed  waveforms  were 
computed  for  all  possible  P  and  T  axis  locations  and  slip  vector  orientations  and  con¬ 
toured  on  a  projection  of  the  focal  sphere.  The  error  contour  at  which  the  synthetic 
waveforms  distinctly  diverged  from  the  observed  waveforms  was  established  by  eye.  The 
procedure  shows  that  the  analysis  of  long-period  waves,  at  least  in  this  case,  provides 
much  better  costraint  on  focal  mechanism  orientation  than  does  first  motion  data  alone. 

In  Section  IV,  results  are  presented  of  modeling  observations  of  S,  SS,  and  SSS 
where  the  event  occurs  near  tectonic  margins  and  are  recorded  om  stable  continents. 
The  modeling  for  laterally  varying  structure  is  essentially  along  a  profile  from  California 
to  Greenland.  The  models  are  allowed  to  be  locally  dipping  with  the  lithosphere  thicken¬ 
ing  with  age  at  the  expense  of  a  dwindling  low  velocity  zone.  Lateral  variation  does  not 
appear  to  be  required  for  depths  greater  than  400  km  along  this  particular  profile.  The 
best  fitting  model  has  a  large  increase  in  lithospheric  thickness  near  the  Rocky  Moun¬ 
tain  Front,  roughly  an  increase  in  75  km  over  a  horizontal  distance  of  400  km  or  less. 
The  low  velocity  zone,  with  a  velocity  of  4.4  km/sec,  is  replaced  by  a  much  faster  upper 
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Modeling  Lg  wave  propagation  in  laterally  heterogeneous  media. 

Janice  Regan 

SUMMARY 

We  have  developed  finite  element  forms  of  the  seismic  Representation  Theorem 
(RT)  to  model  Lg  waves  crossing  continental  margins.  We  have  verified  the  accuracy  of 
this  formulation  using  a  layer-over-halfspace  model  for  which  a  solution  can  be  indepen¬ 
dently  determined  using  conventional  mode  sum  calculations.  We  also  used  this  test 
case  to  illustrate  the  efficiency  of  an  absorbing  boundary  condition  designed  to  reduce 
computation  time.  In  the  modeling  of  transition  regions,  the  decrease  in  crustal  thick¬ 
ness  from  continents  to  oceans  focuses  seismic  energy.  This  results  in  surface  amplitudes 
that  are  50%  larger  in  the  oceanic  crust  near  the  end  of  the  transition  than  those  in  the 
adjacent  continent.  However,  the  amplitude  of  the  Lg  waves  rapidly  attenuates  as  they 
propagate  to  larger  distances  in  the  oceanic  crust.  Conversely,  amplitudes  are  reduced 
when  the  Lg  wavetrain  passes  through  the  transition  from  ocean  to  continent  because 
energy  disperses  into  the  thickening  crustal  layer. 

INTRODUCTION 

The  principal  advantage  of  our  RT  technique  for  propagating  Lg  waves  is  that  it 
enables  us  to  perform  the  calculations  for  mode  sums.  We  are  not  restricted  to  mode  by 
mode  calculations  through  the  transition  region  as  in  previous  studies.  In  the  results 
reported  here,  Lg  is  represented  by  the  fundamental  and  first  5  higher  SH  modes.  These 
6  modes  are  then  simultaneously  propagated  through  the  Finite  Element  (FE)  grid 
representing  the  transition  region  of  interest.  A  further  advantage  coupling  is  that  the 
FE  calculations  are  driven  by  specifying  the  displacement  time  histories  at  a  column  of 
nodes  which  represents  the  boundary  between  the  laterally  homogeneous  region  and  the 
transition  zone.  These  time  histories  can  be  calculated  using  standard  techniques  when 
the  source  is  separated  from  the  margin  by  long  continental  paths.  Alternatively,  they 
can  be  input  from  other  FE  calculations,  as  in  the  case  of  the  continent-ocean-continent 
transition  discussed  below.  We  should  point  out  however  that  our  technique  is  only 
valid  for  waves  propagating  at  normal  incidence. 
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absorbing  were  also  calculated  using  the  short  grid  for  the  FE  portion  of  the  calculation. 
However,  in  this  case  the  absorbing  BC  was  imposed  at  the  dotted  vertical  line  in  Figure 
la.  The  results  show  that  the  absorbing  BC  is  very  efficient  for  Lg  mode  sum  input. 

Reflections  from  the  bottom  edge  of  the  grid  should  also  be  considered.  In  the  case 
of  an  SH  pulse  propagating  in  a  half-space  reflections  from  the  bottom  of  the  grid  also 
introduce  difficulties.  However,  for  the  case  of  Lg  mode  sums  propagating  in  a  layer 
over  a  half-space,  the  Lg  wave  interacts  predominantly  with  the  real  boundary  between 
the  layer  and  the  half-space.  Introducing  absorbing  BC’s  on  the  bottom  of  the  grid 
makes  no  perceptible  difference  in  the  results.  Since  including  the  BC  slightly  increases 
the  number  of  operations  required  at  each  boundary  node  at  each  time  step,  the  bottom 
edge  of  the  grid  is  not  constrained  in  any  of  the  examples  described  here. 

The  FE  grids  used  to  study  the  propagation  of  an  Lg  wave  along  a  continental  to 
oceanic  to  continental  path  are  shown  in  Figure  2a.  The  forward  transition  that  takes 
the  energy  across  a  continental  to  oceanic  transition  region  is  modeled  as  a  smooth  thin¬ 
ning  of  the  crustal  layer  and  a  gradual  introduction  of  a  water  layer  over  a  distance  of 
twenty  kilometers.  The  pair  of  vertical  lines  at  the  left  end  of  the  grid  represent  nodes 
constrained  to  move  with  the  input  Lg  displacement  histories.  The  2  vertical  lines 
labeled  Fcl  and  Fc2  show  receiver  sections  used  to  determine  the  effect  of  the  transition 
on  the  distribution  of  energy  with  depth.  The  pair  of  vertical  lines  labeled  A  are  nodes 
for  which  displacement  histories  are  recorded  to  use  as  input  to  the  reverse  transition. 
The  reverse  transition  grid  takes  the  wave  across  an  oceanic  to  continental  transition 
region.  This  transition  is  modeled  as  a  smooth  thickening  of  crust  and  thinning  of  the 
water  layer  over  a  distance  of  twenty  kilometers.  The  pair  of  vertical  lines  at  the  left 
end  of  the  grid  are  the  nodes  constrained  to  move  with  the  Lg  displacement  histories 
recorded  at  the  column  marked  A  in  the  forward  transition  grid.  The  2  vertical  lines 
marked  Rcl  and  Rc2  show  receiver  sections  used  to  determine  the  effect  of  the  transi¬ 
tion  on  the  energy  distribution  with  depth.  Dots  on  Fcl,  Fc2,  Rcl,  and  Rc2  indicate 
positions  of  receivers  for  which  displacement  histories  are  illustrated  in  Figure  2b.  Open 
circles  refer  to  nodes  for  which  displacement  histories  are  plotted  in  Figure  2c. 

The  passage  of  the  Lg  wave  through  the  transition  region  effects  the  distribution  of 
energy  with  depth.  The  energy  impinging  on  the  crust/water  sloping  boundary  is 
mostly  reflected.  Since  no  energy  is  transmitted  into  the  water  layer,  the  remaining 
energy  is  focussed  along  the  sloping  boundary  into  the  oceanic  crustal  layer.  This  can 
be  seen  from  Figures  lb  and  2b.  Comparing  the  seismograms  in  row  Fcl  of  Figure  2b 
and  the  synthetic  traces  for  a  layer  over  a  half-space  shown  in  Figure  lb,  the  difference 
in  wave  shape  indicates  interference  of  a  reflected  wave  and  the  propagating  Lg. 
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in.  A  Body-Wav*  Analysis  of  tha  1966  Gisborna,  Naw 

Zoaland,  Earthquake. 

Terry  H.  Webb",  S.  G.  Wesnousky  and  Donald  V.  Heimberger. 

Setsmologlcal  Laboratory, 

California  Institute  of  Technology, 

Pasadena,  California  91 125. 


Abstract 

The  Ml  -  6.2  Gisborne  earthquake  of  March  4,  1966  occurred  along  the  East 
Coast  of  the  North  Island,  New  Zealand.  Modeling  of  P  and  S  body  waves  shows  that 
the  focal  mechanism  of  this  event  is  consistent  with  northwestward  thrusting  of  the 
Pacific  Plate  beneath  the  North  Island  ($p =249°, <5  =  25°  and  \=131°).  The  focal 
depth  is  constrained  to  1 8  km,  significantly  less  than  the  values  of  25  to  30  km  com¬ 
puted  from  local  network  data.  Estimates  of  the  scalar  moment,  source  duration  and 
stress-drop  for  the  event  are  4x1024  dyne-cm,  2  to  3  seconds,  and  20  to  120  bars 
respectively.  Cross-correlation  errors  of  synthetic  to  observed  waveforms  were  com¬ 
puted  for  all  possible  P  and  T  axis  locations  and  slip  vector  orientations  and  con¬ 
toured  on  a  projection  of  the  focal  sphere.  The  error  contour  at  which  the  synthetic 
waveforms  distinctly  diverged  from  the  observed  waveforms  was  established  by  eye. 
The  procedure  shows  that  the  analysis  of  long-period  body  waves,  at  least  in  this 
case,  provides  much  better  constraint  on  focal  mechanism  orientation  than  does  first 
motion  data  alone. 


st:-0»uptvy»fci  Otvtrtan,  O.  Sot  13 10,  Wellington, -Ms-w  Julinu. 


3- 


ftecent  reanalysis  of  the  New  Zealand  network  data  does,  however,  show  the  depth 
of  the  earthquake  to  be  welt-constrained  at  25  to  30  km  (E.  G.  C.  Smith,  personal 
communication  1984). 

Teleseismie  Data 

The  WWSSN  stations  that  recorded  P  and  S  waves  for  the  Gisborne  earthquake 
are  shown  in  Figure  2.  The  azimuthal  distribution  of  stations  recording  the  event  is 
poor,  primarily  due  to  the  relatively  small  size  of  the  event  and  the  wide  expanse  of 
the  Pacific  Ocean  which  lies  to  the  north,  east  and  south  of  New  Zealand.  In  summary, 
8  reliable  P  waveforms  and  2  SH  waveforms  were  recovered  from  the  WWSSN  library. 
The  P-waves  are  quite  similar  to  within  the  noise  level  which  reflects  their  close  spa¬ 
tial  distribution  on  the  focal  sphere.  Some  differences,  such  as  the  amplitude  of  the 
clear  second  arrival  and  of  the  final  backswing,  are  not  consistent  with  relative  sta¬ 
tion  location,  so  probably  do  not  reflect  properties  of  the  source.  Interestingly,  none 
of  the  stations  that  recorded  P-waves  showed  evidence  of  SH  energy,  possibly  a 
result  of  being  situated  near  a  S-wave  node.  Although  TAU  (4=23.5°)  and  RIV 
(4=22.0°)  recorded  tangential  motion,  the  waveforms  were  dominated  by  Love 
waves  and  were  consequently  not  useful.  It  should  also  be  noted  that  the  TAU  P- 
wave  has  been  used  without  modeling  the  possible  effects  of  a  triplication.  The  simi¬ 
larity  of  the  P  waveform  recorded  at  TAU  to  P  waveforms  recorded  nearby  on  the 
focal  sphere  (Figure  2)  suggests  that  the  TAU  P  waveform  has  not  been  complicated 
by  triplication  effects.  The  TAU  waveform  Is  actually  the  radial  component,  as  the 
vertical  one  was  not  recorded. 
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main  effect  of  lowering  the  surface  velocities  was  to  reduce  the  size  of  the  sP 
phase. 

Results 

Synthetic  waveforms  for  the  best-fitting  solution  (^=249®  ,5=25°  and  \=131a) 
are  compared  to  the  observed  waveforms  in  Figure  2.  All  of  the  P-wave  data  show  a 
prominent  depth  phase  about  10  seconds  after  the  first  arrival.  Some  difficulty  was 
initially  experienced  in  deciding  whether  this  was  the  phase  pP  or  sP.  As  this  radically 
affects  the  final  best-fitting  mechanism  and  our  result  for  the  focal  depth  differs 
from  previous  studies,  we  examine  this  question  in  detail.  Interpretation  of  the  strong 
depth  phase  as  pP  is  consistent  with  the  25  to  30  km  depth  determined  from  regional 
network  data.  However,  the  best-fitting  mechanism,  assuming  a  28  km  focal  depth 
produces  synthetic  waveforms  that  do  not  fit  well  at  three  stations  (GSC,  PAS  and 
TAU).  Relatively  good  fits  to  the  data  from  all  stations  were  achieved  if  synthetics 
were  made  with  an  1 8  km  focal  depth,  which  implies  that  the  prominent  depth  phase 
is  sP,  not  pP.  The  synthetics  for  each  depth  for  two  of  the  sensitive  stations,  GSC 
and  TAU,  are  shown  in  Figure  3.  The  synthetic  SH  pulse  width  for  an  assumed  focal 
depth  of  28  km  is  wider  than  that  observed  at  GSC,  whereas  a  good  match  is 
obtained  when  an  18  km  depth  is  used.  Since  the  width  of  the  SH  waveform  is 
governed  by  the  S-sS  time  delay,  this  provides  an  unambiguous  constraint  on  the 
depth.  In  addition,  the  choice  of  an  18  km  focal  depth  provides  a  superior  fit  to  the 
TAU  waveform  compared  to  that  for  the  28  km  depth.  At  the  greater  depth  there  is  a 
strong  trade-off  between  the  TAU  P  and  GSC  SH  waves  becoming  nodal. 

Further  evidence  for  identifying  the  prominent  depth  phase  as  sP  was  found  by 
comparing  short-period  records  from  the  South  American  stations  ARE  and  PEL,  for 
which  there  are  no  long-period  data,  with  short-period  records  from  two  stations  at 
opposite  azimuths  -  RAB  and  SNG.  The  vertical  component  waveforms  are  shown  in 
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A  trapeziodal  time  function  of  0.7,  1.0  and  0.7  seconds  was  used  for  the  syn¬ 
thetics.  Other  shapes  couid  fit  equally  well,  but  the  overall  source  duration  was 
determined  to  be  in  the  range  2  to  3  seconds.  This  is  within  the  range  of  source  dura¬ 
tions  (for  an  event  of  this  size)  for  the  group  of  plate  boundary  events  compiled  by 
Liu  and  Kanamori  (1980).  Assuming  a  circular  fault,  rays  leaving  at  45s  to  the  normal 
to  the  fault  plane  and,  using  the  formula  derived  by  Getler  (1976)  and  Ebel  et  al. 
(1978),  we  calculated  a  fault  radius  of  2.6  to  4.0  km.  The  stress-drop,  calculated 
from  Brune's  (1970)  circular  fault  model,  is  then  20  to  1 20  bars. 

The  most  subjective  part  in  determining  the  degree  of  constraint  on  the  final 
focal  mechanism  was  in  assessing  for  which  values  of  cross-correlation  error  the 
synthetic  fits  became  unacceptable.  Once  this  was  established  by  eye,  an  algorithm 
based  on  the  program  FOCPLT  (Whitcomb,  1973),  which  tests  all  possible  double  cou¬ 
ple  mechanisms  in  3.5  degree  steps,  was  used  to  determine  the  cross-correlation 
errors  for  all  mechanisms,  P  and  T  axis  locations  and  slip  vector  orientations.  These 
could  then  be  contoured  to  the  level  of  satisfactory  fits  and  could  also  be  compared 
with  similar  plots  made  using  the  first  motion  data. 

Figure  5  shows  the  range  of  possible  mechanisms  determined  from  body-wave 
modelling  (A)  and  first  motion  data  (B),  while  Figure  6  shows  possible  slip  vector 
orientations  and  P  and  T  axis  locations  corresponding  to  this  range  of  mechanisms  for 
the  event.  The  dashed  curves  in  Figure  6  show  the  possible  range  using  first  motion 
data  while  the  solid  curves  denote  the  range  of  acceptable  values  that  produce 
satisfactory  agreement  of  synthetic  and  observed  waveforms.  The  body-wave 
modeling  appears  to  have  reduced  the  uncertainty  In  slip  vector  and  P  and  T  axis 
locations  by  about  a  factor  of  two,  and  is  also  more  reliable  since  more  than  just  first 
motions  are  being  matched. 

The  failure  to  match  the  final  backswing  at  four  stations  (PMG,  CTA,  HKC  and 
SNG)  with  the  synthetics  could  suggest  some  source  complexity,  such  as  a  second 
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The  possible  range  of  horizontal  projections  of  slip  vectors  for  our  solution 
favour  slip  in  the  azimuth  range  265s  to  335°.  The  relative  plate  motion  in  this  region 
is  directed  at  an  azimuth  of  270s  (Walcott,  1978).  Since  this  is  near  one  extreme  of 
the  range  of  slip  vectors,  It  appears  likely  that  slip  occurred  more  nearly  perpendicu¬ 
lar  to  the  strike  of  the  Wadati-Benioff  zone  than  the  plate  motion  direction  and  that 
the  remaining  component  of  plate  motion  is  being  accommodated  by  strike-slip  move¬ 
ment  in  a  shear  zone  further  inland,  as  proposed  by  Walcott  (1978). 

Conclusions 

Analysis  of  P  and  S  waveforms  has  allowed  us  to  provide  more  constraint  on  the 
focal  mechanism  of  the  1 966  Gisborne  earthquake  than  has  previously  been  afforded 
by  first  motion  data  alone.  Similar  analysis  of  other  moderately  sized  earthquakes 
along  the  New  Zealand  plate  boundary  may  help  to  refine  and  constrain  models  of 
plate  motion  and  deformation. 

The  preferred  strike,  dip,  rake  and  depth  of  the  1 966  Gisborne  earthquake  are 
249°,  25s,  131°  and  18  km,  respectively.  The  simplest  Interpretation  of  this  earth¬ 
quake  is  that  it  is  a  typical  subduction  zone  thrust  event  with  a  small  right-lateral 
component  due  to  the  oblique  convergence  of  the  Pacific  and  Indian  plates  in  this 
region.  There  is  no  evidence  from  our  study  to  suggest  that  this  region  is  currently 
undergoing  extension  normal  to  the  strike  of  the  subducted  plate  as  suggested  by 
strain  data  for  the  period  1931  to  1970  (Walcott,  1978),  although  there  may  be  a 
difference  between  the  strain  occurring  in  the  upper  crust  and  that  occurring  near 
the  plate  interface.  Determination  of  other  reliable  earthquake  mechanisms  can  help 
to  clarify  these  observations  by  showing  whether  the  Gisborne  event  occurred  during 
a  short  period  of  east-west  compression  in  a  region  that  has  apparently  otherwise 
had  steady  east-west  extension  in  the  past  50  years,  or  whether  there  is  a  change 
in  the  direction  of  principal  horizontal  compression  with  depth  in  the  crust.  It  is 
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mechanism  with  <p=2AQa,  6=2Sa,  \=131°  and  H=18  km,  while  the  bottom 
traces  are  for  <p=310°,  6=35°,  X=105®  and  H*28  km. 

Figure  5.  Possible  range  of  mechanisms  for  the  body-wave  solution  (A)  and  for 
short-period  first  motion  data  (8)  on  equal-area  lower  hemisphere  projections. 

Figure  6.  Contours  of  possible  slip  vector  orientations  (A)  and  P  and  T  axis  loca¬ 
tions  (B)  for  the  body  wave  data  (solid  lines)  and  short-period  first  motion 
data  (dashed  lines).  Projections  are  as  in  Figure  5. 
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1.  Introduction 

In  a  recent  paper,  (Grand  and  Helmberger  (1984a)),  upper  mantle  models  appropri¬ 
ate  for  pure-path  tectonic  and  shield  regions  were  presented.  By  analyzing  the  (SS-S) 
travel  times  and  wraveform  information  it  became  apparent  that  the  seismograms  could 
be  partitioned  into  these  two  groups  and  interpreted  in  terms  of  vertical  structure,  see 
Fig.  1.  Model  SNA  was  derived  from  Canadian  shield  observations  but  has  proven  use¬ 
ful  in  modeling  observations  from  other  shields,  Rial,  et  al.  (1984).  The  TNA  model 
appears  appropriate  for  younger  oceanic  structures  and  the  more  tectonic  parts  of 
Western  North  America.  Older  oceans  have  a  thicker  lithosphere  approaching  the  physi¬ 
cal  characteristics  of  the  shield  but  conforming  to  the  TNA  model  with  depth  as 
displayed  in  Fig.  1.  Model  ATL  was  derived  for  the  100  M  year  old  portion  or  the 
Western  Atlantic,  Grand  and  Helmberger  (1984b). 

In  this  paper,  we  will  investigate  the  nature  of  the  horizontal  boundaries  that  con¬ 
nect  the  TNA  and  SNA  provinces  by  modeling  the  seismograms  that  sample  both  areas. 

The  events  used  in  this  study  occur  along  the  western  seaboard  with  recordings 
taken  from  the  WWSSN  and  CSN  stations.  The  locations  of  the  stations  and  events  are 
presented  in  Fig.  2.  The  events  are  also  listed  in  Table  1.  The  reflection  points  at  the 
free  surface  for  the  SS  phases  occur  mostly  in  the  midwest  with  the  bottoming  or  turn¬ 
ing  points  located  in  the  Western  United  States  and  in  Eastern  Canada.  The  lateral 
structure  across  this  region  is  probably  three-dimensional  but  the  general  trend  is  from 
tectonic  to  shield.  Thus,  to  simplify  a  complex  situation,  we  will  assume  that  the  data, 
associated  with  the  station  locations  displayed  in  Fig.  2,  can  be  averaged  from  north  to 
south,  i.e.  observations  at  SHA  of  C2  can  be  compared  to  MDS  of  05,  etc.  In  general, 
these  seismograms  do  share  common  properties  as  pointed  out  earlier  in  Grand  and 
Helmberger  (1984a).  Figure  3  displays  representative  seismograms  from  such  paths  and 
synthetics  appropriate  for  the  models  TNA  and  SNA.  Examining  this  figure  we  see  that 
the  data  are  closer  to  the  SNA  synthetics  in  timing  but  are  more  similar  to  the 
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by  the  two  methods  agree  about  as  well  as  in  the  flat-layered  case,  see  Grand  and  Helm- 
berger  (1984a).  The  first  arrival  in  the  GRT  result,  is  relatively  larger  at  the  nearest  dis¬ 
tances  due  to  tunnelling  through  the  thin  lid.  At  distances  near  28°  the  second  negative 
peak  in  the  WKB  synthetics  is  too  large  compared  to  the  GRT  results.  This  arrival  is 
the  back  branch  of  the  600  km  triplication  and  such  arrivals  tend  to  be  too  large  near 
the  triplication  tip,  see  Burdick  and  Orcutt  (1978)  for  example.  The  last  arrivals 
denoted  by  the  arrows  are  truncation  effects.  Both  can  be  eliminated  or  modified  by 
adding  more  generalized  rays,  see  Helmberger  et  al  (1985).  However,  we  are  primarily 
interested  in  mapping  the  position  of  the  triplications  and  can  tolerate  some  synthetic 
deficiencies.  We  will  also  omit  the  arrivals  guided  by  the  growing  lithosphere  after  giv¬ 
ing  a  brief  review  of  their  properties. 

A  profile  of  synthetics  containing  the  complete  generalized  ray  set  for  the  linear 
model  is  displayed  in  Fig.  6.  As  in  a  recent  paper  by  Grand  and  Helmberger  (1984a)  we 
have  broken  up  the  rays  into  two  groups,  those  that  travel  in  the  crust  and  lid  or  litho¬ 
sphere,  and  those  that  return  from  below  the  low  velocity  zone,  LVZ.  All  the  responses 
are  normalized  to  the  top  trace.  The  short  period  precursors  to  the  long  period  Love 
wave  are  multiples  developing  in  the  growing  lid.  Upper  mantle  arrivals  become  contri¬ 
butors  to  the  complete  response  at  about  18.2°(1800  km)  for  this  particular  model.  The 
two  sharp  spikes  occurring  in  this  trace  are  the  onset  of  the  400  km  triplication,  essen¬ 
tially  S  and  sS.  Rays  bottoming  between  200  and  300  km  produce  the  main  arrivals. 
Synthetic  seismograms  assuming  a  t/  —  3  sec  and  a  far-field  trapezoidal  time  history  of 
(1,1,1)  are  given  on  the  right.  We  probably  over  attenuated  the  lid  arrivals  in  these  syn¬ 
thetics  since  when  we  can  see  these  precursors  on  long  period  observations  they  are  also 
strong  on  the  short  period  records  implying  a  high  Q. 

We  do  not  see  precursors  for  events  near  the  western  portion  of  the  United  States, 
see  Grand  and  Helmberger  (1984a),  although  they  can  be  observed  in  Canada  for  some 
events  located  near  the  edge  of  the  shield.  Since  we  will  be  studying  mostly  western 
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Comparison  with  Observations 

Selecting  a  set  of  observations  to  use  for  comparison  with  the  synthetics  discussed 


A 


in  the  previous  section  may  seem  easy  given  the  large  numbers  of  events  available.  How¬ 
ever,  because  of  radiation  pattern  effects  we  may  not  obtain  simple  SH  pulses  along 
desired  paths.  SV  coupled  PL  waves  make  working  with  rotated  records  difficult  in 
many  situations  as  pointed  out  earlier,  Helmberger  and  Engen  (1974).  Thus,  we  chose 
only  observations  that  are  relatively  naturally  rotated  or  have  stations  located  on  SV- 
nodes.  The  latter  condition  can  be  determined  by  studying  the  event  teleseismically  and 
since  we  use  only  strike-slip  events  the  nodes  for  P  and  SV  occur  along  the  same 
azimuth.  All  the  events  used  in  this  analyses  are  relatively  shallow,  less  than  15  km. 
The  standard  depth  used  in  all  the  synthetics  was  taken  as  8  km,  unless  stated  other¬ 
wise.  A  map  displaying  the  stations  and  events  used  in  the  direct  S  waveshape  analyses 
satisfying  the  above  criterion  is  given  in  Fig.  2  as  discussed  earlier. 

A  comparison  of  some  of  these  observations  with  the  best  fitting  model  displaying 
the  onset  of  the  400  km  triplication  is  given  in  Figs.  13  and  14.  Many  of  these  events 
are  quite  small  and  have  shorter  time  functions  than  the  (1,1,1)  assumed  in  Fig.  11,  the 
event  C2  with  M  =  5.2  being  a  good  example.  The  synthetics  in  these  figures  were  com¬ 
puted  with  the  WKBJ  code  and  thus  do  not  contain  the  tunnelled  energy  which  appears 
as  a  long  period  onset  to  the  mantle  S  phase  as  displayed  in  Fig.  15.  The  effect  disap¬ 
pears  at  larger  ranges  as  the  S  phase  leaves  the  shadow  zone  boundary,  in  this  case  near 
17°.  We  interpret  the  gradual  beginning  of  the  observed  data  at  DAL  and  GOL,  as  pri¬ 
marily  due  to  this  tunneled  energy.  These  same  events  as  recorded  at  the  more  tectonic 
type  stations,  TUC  and  ALQ,  have  sharp  onsets  with  little  evidence  of  any  precursors. 
However,  the  thickness  of  the  lithosphere  is  likely  to  be  highly  variable  across  the 
tectonic-shield  transition  and  given  the  uncertainty  in  attenuation  effects  and  source 
durations  we  do  not  think  that  detailed  modeling  of  the  precursor  is  justified  at  this 
stage.  Nevertheless,  the  400  km  triplication  is  obvious  in  the  data  and  this  feature  alone 
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pattern  occurred  at  ranges  31  to  35°  in  the  pure  tectonic  model.  These  two  arrivals,  in 
Fig.  18,  merge  at  larger  ranges  producing  a  strong  peak.  The  400  km  and  600  km  tripli¬ 
cations  cross  near  45°.  The  ratio  of  SS  to  S  becomes  very  large  here  but  depends  criti¬ 
cally  on  the  source  duration  and  the  precise  positions  of  the  triplications. 

It  appears,  from  this  profile,  that  the  SS  phase  for  our  model  is  relatively  fast  out 
to  about  38°.  This  would  indicate  that  the  lithosphere  is  possibly  too  thick  at  epicentral 
distances  less  than  15°,  and  should  grow  more  slowly  before  then,  with  perhaps  a  sharper 
jump  at  the  Rocky  Mountain  Front.  Other  profiles  of  data  along  this  azimuth  are 
displayed  in  Grand  and  Helmberger  (1984a).  Some  of  these  events  are  located  far 
enough  to  the  South  to  shift  the  interference  pattern  with  respect  to  the  stations. 

There  is  a  tendency  for  the  ratio  of  S  to  SS  amplitudes  to  vary  somewhat  from 
event  to  event.  This  feature  of  the  data  has  not  been  addressed  here  but  probably  is 
related  to  small  amounts  of  dip-slip  components  present  in  these  events.  Corrections  for 
such  effects  would  be  required  if  the  ratios  were  being  used  for  attentuation  measure¬ 
ments. 

The  other  two  models  produce  reasonable  fits  to  this  profile  although  the  synthetics 
for  SS  are  too  slow.  A  comparison  of  two  keys  distances  is  given  in  Fig.  19  where,  as  in 
the  direct  S  case,  the  FF  model  appears  slightly  too  early,  particularly  at  A  =  36.5°.  At 
A  =  47.6°  the  FF  synthetic  fits  the  observed  waveform  quite  well.  Note  that  the  SFA 
observation  is  from  event  C3,  at  the  California-Baja  border  while  the  observation  at 
SCH  is  of  the  Bl  event  near  the  tip  of  Baja.  The  first  leg  of  SS  samples  a  zone  near 
GOL  for  the  former  path  and  near  JCT  for  the  latter  path.  See  figure  2.  From  many 
detailed  comparisons  between  the  various  observations  and  proposed  models,  we  reach 
the  conclusion  that  substantial  three-dimensional  structure  is  present  in  the  data.  How¬ 
ever,  the  RMFS  model  appears  to  be  a  good  average  two-dimensional  model. 

A  cross-section  of  velocity  structure  with  the  SS  ray  paths  appropriate  for  those 
arrivals  bottoming  near  the  400  km  discontinuity  is  given  in  Fig.  20.  Note  that  the  ray 
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Discussion 

The  RMFS  model  fits  a  relatively  large  amount  of  travel  time  and  waveform  data 
and  implies  a  rapid  horizontal  gradient  in  velocity  near  the  Rocky  Mountain  Front. 
Exactly  where  the  model  has  this  sharp  change  relative  to  the  surface  geology  is  uncer¬ 
tain  because  of  the  averaging  process  discussed  earlier.  Nevertheless,  it  is  interesting  to 
compare  the  predicted  delay  time  produced  by  this  model  with  the  observed  travel  time 
residuals,  see  Fig.  22.  We  suppose  that  stations  BKS  and  GSC  are  near  the  western 
edge  of  the  model  and  plot  the  relative  travel  time  delay  appropriate  for  vertical 
incidence.  Correcting  the  data  to  vertical  incidence  would  reduce  the  total  variation  by 
about  15%  so  that  the  over-all  agreement  between  the  model  and  data  is  quite  good. 
Note  that  the  most  rapid  changes  in  the  data  occur  near  stations  RCD  and  DAL  with 
about  a  4  sec  shift  in  travel  times.  These  stations  show  considerable  scatter  with  strong 
evidence  of  azimuthal  dependence.  Some  of  this  scatter  is  probably  caused  by  slab 
extension  in  the  source  region,  Jordan  (1977),  and  anomalies  in  the  lower  mantle,  Lay 
(1983).  However,  the  analysis  presented  here  would  not  eliminate  the  possibility  of  hor¬ 
izontal  gradients  sharper  than  those  suggested  by  RMFS  and  if  such  strong  transitional 
boundaries  exist  we  might  expect  to  see  some  particularly  interesting  travel  time  pat¬ 
terns  along  the  boundary.  Note  that  such  geologic  boundaries  would  probably  not  be 
linear  as  supposed  in  this  study.  A  manifestation  of  such  structure  on  the  character  and 
waveshapes  of  recorded  teleseismic  signals  has  not  been  noted  in  the  literature  but,  on 
the  other  hand,  the  data  and  methods  required  to  synthesize  such  effects  are  only  now 
becoming  available. 

There  is  a  tendency  for  the  S  residuals  for  stations  near  the  Eastern  seaboard  to 
become  less  negative.  For  instance,  the  values  observed  at  Bermuda,  BEC,  appear  to  be 
only  slightly  smaller  than  predicted  by  JB,  a  result  completely  compatible  with  the  ATL 
model  presented  in  Fig.  1,  see  Grand  and  Helmberger  (1984b). 
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investigations  of  remote  regions  where  little  direct  data  is  available. 
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arrivals. 

Figure  16.  Observations  and  synthetics  displaying  the  interaction  with  the  600  km 
triplication.  The  various  time  histories  are  indicated  for  each  synthetic. 

Figure  17.  Reduced  travel  times  from  observations  compared  with  WKBJ  synthetics 
for  possible  models. 

Figure  18.  Comparison  of  (SS-S)  observations  and  synthetics.  Unmarked  seismo¬ 
grams  are  from  event  Bl.  Others  have  event  numbers  as  in  Table  1.  Timing 
lines  for  SS  are  at  same  time  for  both  data  and  synthetics. 

Figure  19.  Comparison  of  two  example  distances  with  corresponding  synthetics. 

Figure  20.  Ray  paths  sampling  the  upper  400  km  velocity  structure  as  a  function  of 
lateral  position. 

Figure  21.  Observations  and  synthetics.  The  dotted  lines  have  been  included  for 
timing  purposes. 

Figure  22.  Comparison  of  relative  JB  residuals  with  synthetics  (RMFS),  modified 
from  Lay  (1983).  The  line  is  dotted  beyond  SCP  since  many  of  these  stations  are 
near  the  Atlantic  coast  and  expected  to  be  late,  see  Grand  and  Helmberger 
(1984b). 

Figure  23.  Possible  velocity  cross-section,  upper  400  km,  from  California  to  Green¬ 
land  on  the  top  with  corresponding  synthetic  on  the  bottom.  The  source  loca¬ 
tions  for  the  two  observations  are  given  below  the  station  names,  see  B4  and  B5 
in  Fig.  2.  The  surface  bounce  points  for  the  SS  phases  are  near  Husdon’s  Bay 
and  South  Dakota,  respectively. 


Figure  4.  Example  of  a  simple  linear  transition  from  tectonic-to-shield  along 
with  ray  paths  appropriate  for  the  two  methods. 


Figure  6.  Step  responses  generated  by  the  crust  and  lid  on  tbe  left,  com¬ 
plete  response  in  the  middle  and  associated  synthetics  on  tbe  right.  A 
trapezoidal  time  function  with  duration  (1,1.1)  was  assumed  for  the  time 
history  and  a  tj  *  3  was  chosen  for  a  guess  at  the  attenuation.  The  dot¬ 
ted  line  indicates  the  mantle  arrival. 
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Figure  22.  Comparison  of  relative  JB  residuals  with  synthetics  (RMFS). 


